Disentangling surface and bulk photoemission using circularly polarized light 
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We show that in the angle resolved photoemission spectroscopy (ARPES) near-surface induced 
fields can be useful for disentangling the surface and bulk related emission. The jump of the 
dielectric function at the interface results in a nonzero term divA in the photoemission matrix 
element. The term happens to be significant approximately within the first unit cell and leads to 
the circular dichroism for the states localized therein. As an example we use ARPES spectra of an 
YBa2Cu307_i crystal to distinguish between the overdoped surface related component and its bulk 
counterparts. 

PACS numbers: 74.25.Jb, 74.72. Hs, 79.60.-i 

Angle resolved photoelectron spectroscopy is subjected 
to a known limitation, which is the surface sensitivity 
caused by a rather small escape depth of emitted photo- 
electrons. According to the universal curve [l| the mean 
free path for the photo electrons in the energy range of 
20-70 eV is about 5-10 A, which is comparable to the 
unit cell of the most of the studied materials. In numer- 
ous photoemission studies it is frequently assumed that 
the spectra from the first unit cell already reflect the bulk 
properties of the material and that the surface related ef- 
fects are of minor influence. However this is not always 
the case. ARPES spectra of the high temperature su- 
perconductor YB2Cu307_5 are just one of the examples 
[3, Q . The signal picked up from the near surface region 
corresponds to an unusually high hole doping level and 
displays no superconductivity, while in the bulk the sam- 
ples are characterized by a narrow superconducting tran- 
sition and uniform doping level, which clearly shows that 
the bulk and surface have different properties. Possible 
way to enhance the bulk contribution consist in the use of 
comparatively small (5-10 eV) 4] or large (^1000 eV)Q 
excitation energies. However, both of these approaches 
have certain restrictions. In the first case it becomes im- 
possible to probe the states in the whole Brilloine zone, 
as /ciimax = V^rnEkhi/h becomes too small, while in the 
high energy case significant deterioration of the momen- 
tum resolution takes place, not to mention that both ap- 
proaches require a specialized light source, which might 
not always be available. 

The fact that the electron escape depth equals A only 
means that the electron intensity is attenuated by the fac- 
tor exp (— z/A), i.e. there are still photoelectrons leaving 
the solid from the depth z > A, but their intensity is 
decreasing according to the exponential law. Therefore 
the spectra do contain a signal reflecting the bulk prop- 
erties, but the problem is how to extract it from under 
the bright surface contribution. Fortunately the bulk and 
surface photoemission differ also due to the pattern of the 
electromagnetic field that excites the electrons and this 
can be used as a "marker" that makes these two spectral 



components discernable. Already in studies of Cu sur- 
face states it was shown that the photoemission matrix 
elements are strongly modified by the term divA, which 
becomes important at the near surface region due to the 
mismatch in the dielectric constant between the solid and 
the vacuum, and the calculated dependence of the matrix 
element on the incidence angle of the linearly polarized 
exciting radiation was found to be in a good agreement 
with experimental data fd\ . The depth down to which the 
term divA modifies the matrix elements depends on the 
field pattern near the surface, exact calculation of which 
turns out to be a very difficult task. Relatively simple 
jellium model employed to describe a free electron metal 
[3] shows that the vector potential of the electromagnetic 
wave A(r) experiences quickly decaying Friedel-like oscil- 
lations with a characteristic scale of a few angstroms. In 
a more recent study of TiS2 , which is closer to the case we 
are going to concentrate later on, the layered structure of 
the crystal was taken into account Q . From these calcu- 
lations it follows that the characteristic thickness dsurf of 
the surface layer is about one unit cell along the normal 
to the surface, since in deeper regions A(r) practically 
stops oscillating reaching its bulk value. 

In this manuscript we show that the discussed term 
divA(r) leads to a circular dichroism for the photoemis- 
sion from the surface layer and how the dichroism can be 
used to obtain information about the states in the bulk 
and at the surface. 

In general, the photoemission matrix elements depend 
on experimental geometry (i.e., mutual position of the 
sample, polarization, wave vector of incident light q and 
the direction of emitted photoelectron k). In Fig 1. we 
depict the geometry of our experiment. The entrance 
slit of the energy analyzer is vertical and its position is 
fixed in space, so that the ARPES spectrum, the so called 
energy-momentum distribution, is just a 2D distribution 
of photo-intensity as a function of energy and the angle rj 
at which electrons enter the analyzer entrance slit. For a 
fixed sample position the wave vector k of the photoelec- 
tron and its projections and ky on the sample surface 
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FIG. 1: (color online). Experimental geometry, (a) linear dependence of the circular dichroism in photoelectron intensity 
D = P'^ — J'^' on the angle rj. The red-white-blue color scale denotes the strength of the dichroism. (b) decomposition of a 
circularly polarized wave in two linearly polarized ones with polarization vectors As and Ac with a phase shift of 7r/4. 



are uniquely defined by the angle rj and can be easily es- 
timated from the drawing. The incident beam Hes in the 
horizontal plane MON with AM ON = /3 = 45°. Differ- 
ent parts of reciprocal space can be probed via rotation 
of the sample around the pivot point O. For instance in 
Fig. lb we show the sample (the crystal primary axes) af- 
ter rotation around the a;-axis by the angle (p, when the 
projection ky becomes negative in contrast to Fig. la, 
where </? = and ky = 0. 

The photoemission matrix element can be estimated 
as a probability of a transition between the initial state 
\i) and the final state |/), which is given by the Fermi 
golden rule [To| : 



(1) 



where the perturbation to the system Hamiltonian 
//port = -^(AV + idivA). The incident circularly 
polarized wave A(r) can be represented as a sum of two 
plane waves (Fig. 1 b): 

A(r) = Ac cos{ujt — qr) -I- As sin(wt — qr). (2) 

It is easy to show that for such a periodic perturbation 



Vc,s = 
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equation can be rewritten as: 



^vacuum ^sample ^vacuum ^^samplc 



Microscopically, as has previously been pointed out, the 
jump in the vector potential appears as a smooth tran- 
sition with a spatial extension c?surf of about one lat- 
tice parameter c, within which the vector potential A(r) 
changes from its vacuum value to the bulk one. For the 
discussed geometry divAs = due to the boundary con- 
ditions, and divAc can be approximated as: 



A sample A vacuum 

divAc « ^ 



(l/e-l)A^ 



CA 



-^surr 
vacuum 



(6) 



C|Ac|sin(/3-(p), 



where C is a complex constant that effectively accounts 
for the width of the transition layer and sample perme- 
ability. To estimate the matrix element (3) we assume 
that the final state is a plane wave: |/) = le*"^""), with k 
being the quasi-momentum of the exited photoelectron. 
Leaving out the irrelevant coefficients yields: 



(i|AcV + i divAc - iAgVIe^'^"-) 



Acik-I- Ask 



1 



■ divA, 



(7) 



At this stage the particular pattern of the near sur- 
face induced field comes into play. Macroscopically, from 
the continuity condition for the components of D and E 
vectors one has: 



E 



vacuum -^sample 



E^ 



vacuum -q vacuum -p^samplc ^-jpsamplc 



: eE , 



(4) 



here subscripts ||/_L indicate vector components paral- 
lel/perpendicular to the sample surface, and e is the sam- 
ple permeability. Using the gauge with the scalar poten- 
tial of the electromagnetic field taken to be zero, the last 



To get the matrix element for the opposite circular po- 
larization one just needs to reverse the direction of the 
vector As, therefore for the dichroism we obtain: 



Anik + Ask 



1 



■ divA, 



Acik — Ask 



1 



■ divA, 



(8) 



Since |Ac«k — Ask| = |Acik + Ask| there will be no 
circular dichroism if divAc equals zero. It is worth to 
mention that this is a rather general statement true for 
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FIG. 2: (color online), (a, b) panels from left to right: position of the energy-momentum cuts in the reciprocal space; sum 
of the spectra measured with opposite circular polarizations; their difference; extracted bulk component; extracted surface 
component. Color scale bars show the absolute intensity. In the block (a) the process of disentangling of the surface related 
anibonding band and the chain band is demonstrated (hv = 50 eV). Block (b) depicts the same procedure, but for the surface 
and bulk components of the bonding band {hv — 55 eV). The antibonding band in this case is suppressed by the unfavorable 
matrix elements, (c) Structure of YBa2Cu307_i near the cleavage plane. 



any geometry as long as the final states can be well ap- 
proximated by the plane waves Q . After simple compu- 
tations, taking into account the experimental geometry, 
the expression (8) reduces to the final form that we are 
going to use: 

I?^ |(*|e'''"-)|'Re(C)|Ac||As||k|sin(77)sin(/3-^). (9) 

An important consequence of this formula is that in the 
described geometry the photoemission signal arising from 
the near surface region would exhibit circular dichroism 
proportional to sin(r/) w rj. It is this particularity that 
can be used to distinguish the contribution to the spec- 
trum arising from the near surface layer, where the term 
with divA is important, from the one coming from the 
deeper regions of the sample. Representing 2D ARPES 
spectra as a sum of a "bulk" and the "surface" compo- 
nent we can write: 

r\r,, Lo) = I^^^'\r,, to) + c.)(l + ar,), 

/"(r;, c^) = c.) -f /^-f(77, c.)(l - ar;), 

where a accounts for the dichroism strength. The con- 
stituent bulk and surface related components can easily 
be obtained: 

/-f(r;,c.)=-^(r'(r;,u;)-/"(,7,^)), 

2a7j ^^^^ 

Now we proceed to the practical application of the de- 
rived formulae using as an example ARPES spectra of 
YBa2Cu307_5 (Fig- 2). The spectra were measured from 
a freshly cleaved surfaces at T = 30 K with SESlOO SCI- 
ENTA electron energy analyzer at the SIS beam line at 



Paul Scherrer Institute. The energy and angular resolu- 
tion were 15 meV and 0.2°, respectively. The excitation 
energy is given in the caption to the figure, further ex- 
perimental details can be found elsewhere 

Similar to many other layered compounds its weak 
kz dispersion allows for a relatively simple exposition 
of ARPES data. In this case the spectral function de- 
pends only on three parameters, so its argument can 
be thought of as a point in the 3D space spanned over 
one energy and two momentum axes. Therefore the 2D 
energy-momentum intensity distributions, already men- 
tioned during the discussion of experimental geometry, 
are practically almost plane cuts through this 3D space 
with the bright features therein corresponding to the 
traces of the renormalized electronic bands [Tlj . In the 
leftmost panels of Fig. 2a, b we show a model of the 
YBa2Cu307_5 Fermi surface map that consists of pairs 
of square-like contours around the S points, correspond- 
ing to the bonding and antibonding bands, and features 
parallel to axis, which are the Fermi level crossings 
of the chain derived band. The red segments denote the 
position in the fc-space of the momentum-energy distri- 
bution given in the next panels, which are: the sum of 
spectra obtained with the opposite circular polarizations; 
the difference of the spectra; the extracted bulk and sur- 
face contributions based on formulae (11). To achieve 
maximal subtraction of the surface component, parame- 
ter a was increased until negative values appeared. 

From the Fig. 2a it follows that the states of the intense 
overdoped antibonding band are supposed to be localized 
at the near surface region, while the chain band has to be 
of bulk origin as it displays no dichroism. To understand 
this one needs to look in detail at the structure of cleaved 
YBa2Cu307_5 shown in Fig. 2c. According to the tun- 
neling experiments [l^, [l3| this crystal cleaves between 
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BaO and CuO layers, which is schematically shown in the 
figure as broken bonds and missing atoms. The remnants 
of the CuO chains, heavily disrupted by the cleavage, are 
unlikely to result in photoemission signal that reminds 
the bulk dispersion of this band. Therefore the signal 
from the nearest to the surface Cu02 bilayer should be 
the brightest. In a view of missing chain structure on top 
of this bilayer, the change of its hole doping level is not 
surprising [3, Q- The absence of the dichroism for the 
chain states simply means that the nearest to the surface 
chains are already out of the region where A(r) strongly 
oscillates resulting in perceivable divA. This also means 
that the next Cu02 bilayer, which we expect to be su- 
perconducting as it is surrounded by the CuO chains on 
its both sides, resides in the region, where divA becomes 
negligible and should exhibit no dichroism, similar to the 
chain band. Indeed, splitting the spectrum of the bond- 
ing band (Fig. 2b) into surface and bulk contributions we 
see that the surface component reminds the spectrum of 
the normal state taken above Tc, exhibiting no unusual 
renormalization. The spectrum of the bulk component 
looks qualitatively different. The strong band renormal- 
ization, which is a known signature of the superconduct- 
ing state [H, [H, [13, d, is clearly visible in the 
spectrum, supporting the expectation. 

Aiming at true quantitative analysis of the disentan- 
gled spectra one needs to be cautious as the dichro- 
ism in photoemission experiments is a rather ubiqui- 
tous phenomenon. Unlike the classical linear magnetic 
dichroism, which rarely exceeds 1% [20], the strength 



(/cr _ /ci)/(/cr _^ jci) ^j^g dcscribcd surfacc related 
dichroism exceeds 60%, so the former can safely be ne- 
glected. However, notable dichroism can also arise from 
the bulk states [21[. While the surface related dichroism 
vanishes at the normal light incidence (Eq. 9, (3 — ip), 
the bulk one is expected to be still non-vanishing if the 
experimental setup (including the sample) possesses def- 
inite handedness and the final states substantially dif- 
fer from the plane waves. To estimate its strength in 
cuprates we can refer to our previous results [l^], where 
the total dichroism for the normal light incidence was 
shown to be not more than 6%. For the arbitrary light 
incidence (23j (when the surface dichroism is vanishing 
due to condition i] = 0) the total dichroism amounts to 
^18%, but for the sample orientation discussed here, i.e. 
when one of the crystal primary axes lies in the MON 
plane, its value is zero again, so that the surface related 
dichroism turns out to be a dominating one. 

To conclude, we have stressed the importance of the 
near surface induced electromagnetic fields for the in- 
terpretation of photoemission data, and pointed out that 
those might be a beneficial factor allowing for distinguish- 
ing between the surface and bulk photoemission. In some 
cases it is also possible to localize the position of the sur- 
face and bulk states more precisely, like in the case of 
YBa2Cu307_5, where the nearest to the surface Cu02 
bilayer happens to be overdoped, while the next bilayer 
can already be treated as the bulk one. The obtained 
results can easily be extended to any other experimental 
geometry. 
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